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Abstract 

We compute the structure of an azamacrocycle, C2qHisNq. Two approx- 
imatively planar elliptical structures with C2 or Cj symmetry are found to 
be nearly degenerate. The roughly circular conformation observed in metal 
complexes turns out to be ~ 0.6 eV higher in energy. We suggest that this 
difference is mainly due to electrostatic interactions. We discuss the results 
on various levels of theory (Hartree-Fock, local density and gradient corrected 
density functional calculations). 
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I. INTRODUCTION 



In recent years, the formation of metal complexes of macrocyclic ligands has become an 
important research field which is rapidly developing. The interest in this field originates from 
the possibility to build supramolecular architectures and thus to finally arrive at tailor-made 
materials with novel properties. The goal to design new materials with a combination of 
desired properties, however, requires a quantitative understanding of the properties of the 
individual units, their response to chemical substitutions as well as their interactions. Of 
prime importance is to understand the molecular structure and the factors which determine 
it. The azamacrocycle, C2&H\%Nq, has been widely used as a ligand, for the formation of 
complexes with a wide variety of metals such as alkali metals!'!, lanthanides!'!, earth alkaline 
met a Is! , lea J, cadmiumi'i or yttrium!. In all these metal complexes, the azamacrocycle 
is approximatively ring-like, with the metal ion bonding to the six nitrogen atoms. It 
was surprising, when the free ligand was synthesized and was found to be approximatively 
planar-elliptical!. 

A full explanation of the difference of the two structures has not yet been given. However, 
modern codes and computers make ab-initio calculations on the Hartree-Fock (HF) and 
density-functional level for these systems possible. We therefore investigated the different 
structures and their electronic properties. The main target of this article is a quantitative 
explanation of the relative stabilities of these structures. Thus we optimized the geometry 
and computed the total energy, charge density and the vibrational spectrum. In addition, 
we try to assess the overall validity of the methods and the approximations involved. We 
compare the results for Hartree-Fock and various density functional calculations to estimate 
the ambiguities introduced by the simplified description of electronic correlations. In the 
following section (0) we give details about the method. The results are presented in section 



III, and finally we summarize the article in section IV. 
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II. METHOD 

All the calculations were done with the code Gaussian 98@. We employ the Hartree- 
Fock method, the local density approximation (LDA) (with the correlation functional V 
from reference [|) and the gradient corrected Hybrid functional B3LYP0. The bulk of the 
calculations was done with a 6-31G Gaussian basis set for C, H and N atoms. Additional 
calculations with an enhanced 6-31G** basis set (i.e. with an additional p-polarization 
function for H and an additional d-polarization function for C and N) were performed to 
investigate the dependence of the results on the choice of the basis set. As an optimization of 
the system in its crystal structure (monoclinic, space group P2i/C) is presently prohibitive, 
we performed a calculation on the isolated molecule. The effects of this approximation will 



be discussed in section |T|. As a first step, a full geometry optimization with respect to 
all parameters and without symmetry was performed. Depending on the geometry at the 
beginning of the optimization process, the final optimized geometry turned out to be either 
C2, Ci or C2V symmetric, the three minima considered in this article. To ensure the stability 
of the optimized structures, in addition vibrational frequencies were calculated. 

III. RESULTS 

The results of the calculations are summarized in tables |,0, [TTT] and [V]. First, we 
discuss the optimized geometries. The first calculated structure is visualized in figure |] using 
Molden§. This structure is approximatively planar-elliptical and has C2 symmetry. Thus, 
it is the one which has the same symmetry as a single C2§Hi%Nq unit of the experimentally 
observed structured (i.e. the inner hydrogens, 41 and 42 in figure p] are on the same face of 
the molecule). The second structure (structure 2) is also approximatively planar-elliptical, 
but has Ci symmetry. In the view in figures 0, it looks virtually identical as structure 1, 
and the difference compared to structure 1 can better be seen in figure |3|. The alternative to 
these flat structures is the approximatively circular structure (structure 3) which is found to 



3 



be another local minimum. We display the latter structure in figure |2j. This latter structure 
is saddle-like; in a view from the top it appears approximatively circular (figure |2|), and, 
when viewing from the side, we note that there is a significant bending (figure [|). 

The energy differences of the structures on the HF, LDA and B3LYP level of theory 
are displayed in table [[[ We find that the saddle-like structure 3 is highest in energy at all 
levels of theory. While HF gives the highest energy difference, this is reduced by LDA and 
even more by B3LYP. The calculation with the best basis set used in this article (6-31 G**) 
gives an energy difference of 0.6 eV at the B3LYP level. Including zero-point vibrations 
hardly changes the energy difference, as the zero-point energy was found to be 0.384 Eh (1 
E h = 27.2114 eV) for structure 1 and 2 and 0.383 E h for the saddle-like structure 3 (B3LYP, 
6-31G** basis set). The computed vibrational frequencies were in the range from ~20 to ~ 
3200 cm" 1 . The energy difference is also consistent with earlier results based on molecular 
mechanics calculationsi, where, depending on the approximations, an energy difference in 
the region from ~ to 0.7 eV was obtained. Comparing the results for structure 1 and 2, 
we note that at the HF level, structure 1 is found to be slightly lower than structure 2. This 
changes when the calculations are performed at the density functional level: now, structure 
2 is consistently lower than structure 1. This energy splitting ranges from +0.0003 Eh to 
-0.002 Eh, The quality of the basis set has a minor influence, and zero-point vibrations are 
negligible for this splitting. 

Although B3LYP calculations are usually reliable, this contradiction of the various levels 
of theory makes the correct prediction of the molecular conformation very difficult. Con- 
cerning the molecular crystal, lattice effects may well be of the order of the magnitude of 
the splitting between structure 1 and structure 2. A safer prediction of the structure of 
the molecular crystal would thus require an optimization of the periodic system which is 
presently prohibitive. 

Structural parameters, together with Mulliken populations are displayed in tables [II], [TTT] 
and [TV]. The charge distributions indicate that all the hydrogen atoms carry charges of ~ 
0.1 |e|, and the nitrogen atoms ~ — 0.5|e|. The carbon atoms carry charges from ~ -0.1 to 
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~ 0.25|e|; obviously they are stronger negative for those atoms with a hydrogen neighbor 
and positive for those atoms with a nitrogen neighbor. The charge distribution is nearly 
consistent for all the structures. 

We further note that the structural parameters vary only weakly when using basis sets 
of different quality. Similarly, the charges varied slightly by about 0.1 |e|. This deviation 
may be viewed as an error bar and the conclusions are not affected. 

A comparison with the experimental structure data was made in table [V|. We used the 
standard deviation of the distance matrix 

/ /i computed jexp\o 

I 2 y^n \ d i,j ~ d i,j ) 

y n(n-l) ^i=l ^3=1 (ZJ? 

as a measure for the error of the computed geometry (with n the number all atoms and 
^computed an( j ^exp ^ e ^ n g com p U t e d and experimental distances between the atoms i and 
j, respectively). We note that the standard deviation is very small both for the structure 
with Ci symmetry as well as the structure with Cj symmetry. However, the structure with 
the C*2 symmetry found experimentally has a slightly higher standard deviation. Finally, 
the structural parameters of structure 3 with Ci v symmetry deviate strongly. 

Comparing the structures, we can thus conclude that electrostatic effects favor the ap- 
proximatively planar structures 1 and 2 over structure 3: the main effect seems to be that 
the two inner hydrogen atoms (41 and 42 in figure [l]) have lowered the energy of the molecule 
by inwards bending and pointing towards the negatively charge nitrogen atoms. This ad- 
ditional attraction would, for example, for the hydrogen 42 and nitrogen 6, be of the order 
0.17 x 0.45 /4 a.u. ~ 0.5 eV and thus seems to be the likely explanation for the preference 
of structure 1. In addition, we note that the hydrogen atoms (33, 50 and 34, 49) in structure 
1 would feel a relatively strong repulsion in a strictly planar structure, and a slight distor- 
tion of the molecule thus lowers the energy. This is indeed found in the experiment and 
in our simulations. A similar distortion is also found in structure 2. Finally, in structure 
3 the hydrogen atoms (39, 43; 41, 45; 40, 44 and 42,46) would feel a repulsion in a planar 
structure; and this repulsion can be reduced by a distortion to a saddle-like structure. 

In summary, we feel that the preference of the planar elliptical structures over the saddle- 
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like can be explained by a purely electrostatic argument. Embedding a positively charged 
metal ion was experimentally found to change the situation: the saddle-like structure is now 
the preferred one. We feel that this effect may have a similar reason: in the approximatively 
planar elliptical structures 1 or 2, the additional repulsion from the two inwards bending 
hydrogen atoms would make it difficult to embed a cation, and moreover, the saddle-like 
structure offers more space to embed a cation. In the saddle-like structure, the metal ion 
feels the electrostatic attraction due to all six nitrogen atoms. However, a full quantitative 
analysis of a system with metal ions will be more difficult because the individual units are 
charged and thus also effects from the anions in the molecular crystal must be included. 



IV. SUMMARY 

We have demonstrated that ab-initio Hartree-Fock and density-functional calculations 
can explain the relative stabilities of flat versus saddle-like conformations of the azamacro- 
cycle C2gHisNq. The charge distribution indicates that electrostatic effects are the reason 
for the energy difference of ~ 0.6 eV between these conformations. However, two flat struc- 
tures with C<i symmetry and Ci symmetry were found to be nearly degenerate. The distance 
matrices of the structures with Ci symmetry and with Cj symmetry are both in very good 
agreement with the distance matrix from experimental data. Finally, we showed that the 
energetically highest structure has a geometry similar to the one observed in systems with a 
metal embedded in the azamacrocycle. We feel that this saddle-like structure is preferable 
for the embedding of a metal ion because of the charge distribution and the available space. 

The main limitation of these calculations was the restriction to molecules instead of 
molecular crystals. New code developments (for example the CRYSTAL code with analytical 
gradients^!) will hopefully make calculations on the molecular crystal feasible, as well as more 
demanding simulations on crystals with embedded metal ions@~i. 
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TABLES 

TABLE I. Relative energies, computed with respect to the energy of the structure with C2 
symmetry, at different levels of theory (1 £^=27.2114 eV) 



method basis set relative energy with respect to 

the structure with C2 symmetry [Eh] 







structure with 


structure with 






Si 1 

Cj symmetry 


Lj2v symmetry 


HF 


6-3 1G 


0.0003 


0.043 


HF 


6-31G** 


0.0009 


0.026 


LDA 


6-3 1G 


-0.0016 


0.038 


LDA 


6-31G** 


-0.0021 


0.028 


B3LYP 


6-3 1G 


-0.0008 


0.033 


B3LYP, with zero-point vibrations 


6-3 1G 


-0.0008 


0.031 


B3LYP 


6-31G** 


-0.0016 


0.022 


B3LYP, with zero-point vibrations 


6-31G** 


-0.0015 


0.022 
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TABLE II. Optimized geometrical parameters and Mulliken populations, B3LYP, 6-31G** 
basis set, approximatively elliptical structure (1) with C2 symmetry (i.e. simultaneously x — > -x 
and y^ -y is a symmetry operation. Only positions of symmetry unique atoms are given. 



n number 


type 


X 






[A] 


1,4 


N 


0.352 


2,5 


N 


-3.187 


3,6 


N 


-2.939 


7,8 


C 


1.376 


9,10 


C 


1.184 


11,12 


c 


-0.122 


13,14 


c 


-1.184 


15,16 


c 


-0.900 


17,18 


c 


-1.948 


19,20 


c 


-4.252 


21,22 


c 


-5.522 


23,24 


c 


-6.692 


25,26 


c 


-6.619 


27,28 


c 


-5.380 


29,30 


c 


-4.190 


31,32 


c 


-2.738 


33,34 


H 


-5.315 


35,36 


H 


-7.521 


37,38 


H 


-7.653 


39,40 


H 


-5.549 


41,42 


H 


-1.569 


43,44 


H 


-2.218 



y z charge 

[A] [A] |e| 

-2.027 -0.144 -0.52 

-1.795 -0.022 -0.53 

1.147 0.275 -0.45 

-2.890 -0.236 0.26 

-4.269 -0.393 -0.11 

-4.765 -0.407 -0.06 

-3.877 -0.295 -0.09 

-2.504 -0.181 0.24 

-1.465 -0.114 0.12 

-0.889 0.112 0.24 

-1.494 0.104 -0.10 

-0.765 0.288 -0.09 

0.613 0.502 -0.10 

1.247 0.498 -0.11 

0.530 0.278 0.27 

2.313 -0.206 0.09 

2.312 0.696 0.09 

1.191 0.680 0.09 

-1.271 0.288 0.09 

-2.569 -0.035 0.10 

-0.443 -0.128 0.18 

-4.201 -0.307 0.12 
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45,46 H -0.302 -5.831 -0.511 0.10 

47,48 H 2.038 -4.935 -0.490 0.10 

49,50 H 3.539 -2.923 -0.650 0.09 
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TABLE III. Optimized geometrical parameters and Mulliken populations, B3LYP, 6-31G** 
basis set, approximatively elliptical structure (1) with Cj symmetry (i.e. simultaneously x — > -x, 
y^ -y and z — > — z is a symmetry operation. Only positions of symmetry unique atoms are given. 



n number 


type 


x 


y 


z 


charge 






[A] 


r t i 

[A] 


r t i 

[A] 


e 


1,4 


N 


0.309 


-1.984 


-0.414 


-0.52 


2,5 


N 


-3.220 


-1.804 


-0.122 


-0.52 


3,6 


N 


-2.926 


1.035 


0.547 


-0.45 


7,8 


C 


1.343 


-2.838 


-0.326 


0.26 


9,10 


C 


1.166 


-4.225 


-0.214 


-0.11 


11,12 


c 


-0.132 


-4.738 


-0.200 


-0.06 


13,14 


c 


-1.204 


-3.857 


-0.261 


-0.08 


15,16 


c 


-0.934 


-2.479 


-0.360 


0.23 


17,18 


c 


-2.005 


-1.461 


-0.358 


0.15 


19,20 


c 


-4.283 


-0.914 


0.055 


0.24 


21,22 


c 


-5.564 


-1.483 


-0.073 


-0.10 


23,24 


c 


-6.721 


-0.752 


0.162 


-0.09 


25,26 


c 


-6.625 


0.580 


0.576 


-0.10 


27,28 


c 


-5.375 


1.172 


0.712 


-0.10 


29,30 


c 


-4.192 


0.464 


0.425 


0.26 


31,32 


c 


-2.702 


2.263 


0.269 


0.09 


33,34 


H 


-5.295 


2.194 


1.070 


0.09 


35,36 


H 


-7.520 


1.151 


0.805 


0.09 


37,38 


H 


-7.693 


-1.224 


0.056 


0.09 


39,40 


H 


-5.611 


-2.531 


-0.350 


0.09 


41,42 


H 


-1.657 


-0.439 


-0.508 


0.15 


43,44 


H 


-2.235 


-4.189 


-0.220 


0.12 
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45,46 H -0.299 -5.808 -0.125 0.10 

47,48 H 2.028 -4.882 -0.139 0.09 

49,50 H 3.490 -2.942 0.088 0.09 
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TABLE IV. Optimized geometrical parameters and Mulliken populations, B3LYP, 6-31G** 
basis set, C2 V symmetric (i.e. x — > -x and y — > -y are symmetry operators). Only positions of 
symmetry unique atoms are given. 



atom number 


type 


X 


y 


z 


charge 


1,2,3,4 


C 


-1.150 


3.105 


0.695 


0.26 


5,6,7,8 


C 


-1.198 


4.323 


1.398 


-0.11 


9,10 


C 


0.000 


4.942 


1.744 


-0.06 


11,12 


N 


0.000 


2.518 


0.343 


-0.47 


13,14,15,16 


C 


2.416 


2.412 


0.359 


0.09 


17,18,19,20 


N 


2.436 


1.347 


-0.342 


-0.39 


21,22,23,24 


C 


3.628 


0.711 


-0.685 


0.22 


or or ot oo 

25,26,27,28 


L 


4.789 


1.395 


-1.079 


n 1 A 

-0.10 


29,30,31,32 


C 


5.930 


0.698 


-1.473 


-0.09 


33,34,35,36 


H 


2.153 


4.764 


1.665 


0.09 


37,38 


H 


0.000 


5.885 


2.283 


0.10 


39,40,41,42 


H 


3.329 


2.863 


0.780 


0.07 


43,44,45,46 


H 


4.771 


2.480 


-1.116 


0.09 


47,48,49,50 


H 


6.810 


1.244 


-1.799 


0.08 
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TABLE V. Error of the distance matrix of the three conformations with respect to the exper- 
imental geometry. 



structure 


method 


basis set 


standard deviation 


(C2 symmetry) 


HF 


6-31G 


0.025 


1 


HF 


6-31G** 


0.029 


1 


LDA 


6-3 1G 


0.029 


1 


LDA 


6-31G** 


0.026 


1 


B3LYP 


6-3 1G 


0.028 


1 


B3LYP 


6-31G** 


0.024 


2 (Q) 


HF 


6-31G** 


0.019 


2 


LDA 


6-31G** 


0.026 


2 


B3LYP 


6-31G** 


0.022 


3 (C 2v ) 


B3LYP 


6-31G** 


0.15 
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FIGURES 

FIG. 1. Optimized planar-elliptical structure with C2 symmetry (structure 1). 
Ci symmetry looks virtually identical in this view. 



Structure 2 with 
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FIG. 2. Optimized saddle- like structure with Civ symmetry (structure 3) 
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